On the determination of CP-even and CV-odd components of a mixed CV Higgs boson 

at e^e" linear colliders 
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We present a method to investigate tiie CV quantum numbers of tiie Higgs boson in tlie process 
e^e~ — > Z<f> at a future e^e~ linear collider (LC), where 4>, a generic Higgs boson, is a mixture of 
CP-even and CP-odd states. The procedure consists of a comparison of the data with predictions 
obtained from Monte Carlo simulations corresponding to the productions of scalar and pseudoescalar 
Higgs and the interference term which constitutes a distinctive signal of CV violation. We present 
estimates of the sensitivity of the method from Monte Carlo studies using hypothetical data samples 
with a full LC detector simulation taking into account the background signals. 



INTRODUCTION 



The future linear e+ e~ collider TESLA is planned to 
work with a maximum center-of-mass energy of 500 GeV, 
extendable to 800 GcV without modifying the original de- 
sign []|. It will have a luminosity of 3.4 x lO'^'* cm~^ s~^, 
a thousand times greater than the LEP at CERN, and 
so it will be well suited for a discovery of a light Higgs 
boson. Even if the Higgs is discovered before at Teva- 
tron (Fermilab) ^ or at the future LHC (CERN, Geneva), 
e + e— colliders are the ideal machines to investigate the 
Higgs sector in the intermediate mass range since all ma- 
jor decay modes can be explored, with the Higgs particle 
produced through several mechanisms [3|- For a light or 
intermediate mass Higgs boson, the Higgstrahlung pro- 
cess e+e" -^ Z<i>, where $ denotes a generic Higgs boson, 
is expected to be the most promising process to study 
its properties and interactions and to search for devia- 
tions from the Standard Model (SM) predictions(see |j| 
and references therein). A comprehensive review of the 
Higgs boson properties has been given in ref. |5j- The 
theory of the Higgs bosons, with emphasis on the Higgs 
scalars of the SM and its non-supersymmetric and su- 
persymmetric extensions has been recently presented in 
ref.|6!|. The spin, parity and charge conjugation quantum 
numbers, J'^'^, of the Higgs boson can potentially be de- 
termined independently of the model. It has been shown 
that measurements of the threshold dependence of the 
Higgsstrahlung cross-section constrains the possibles val- 
ues J^^ of the state |3| • In the minimal Standard Model 
the Higgs mechanism requires only one Higgs doublet 
to generate masses for fermions and gauge bosons 8j. 
It leads to the appearence of a neutral CT^-even Higgs 
{H). In the two-doublet Higgs model (2DHM) or the 
supersymmetric extension of the SM |3|, neglecting CV 
violation, there are two CP-even states (/i, H) and one 
CT'-odd state {A), plus a pair of charged Higgs bosons 



(H^). In a general 2DHM the three neutral Higgs bosons 
could correspond to arbitrary mixtures of CV states and 
their production and decay exhibits CV violation. The 
angular distributions of the Higgsstrahlung cross section 
depends upon whether the $ is CP-even, CP-odd, or a 
mixture [J, |l^ |lj [ij, [iJl • Also the angular distribution 
of the fermions in the Z ^ ff from Z^ production re- 
flects the CV nature of the state * H Efl El El- An 
analysis of the angular distributions of the final state 
fermions in the Higgsstrahlung process with the formal- 
ism of optimal variables has been performed in ref. |l5| . A 
fit to double-differential angular distribution in the pro- 
duction and decay angles results in a clean separation 
between a scalar and pseudoscalar states assuming that 
the Z^ cross section is independent of the CV nature of 
the [13 • Recently, the prospects for the meaurement 
of the pseudoscalar admixture in the hrr coupling to a 
SM Higgs boson was presented |17| . 

In this paper, we present an alternative method that 
simultaneouly uses the distributions of the production 
and decay angles to distinguish the SM-like Higgs boson 
from a CP-odd *" state A, or a CP-violating mixture 
<&. We perform an analyisis of Monte Carlo events that 
takes into account the signals and background, as well as 
a simulation of a TESLA detector response. In the next 
section we shall present the theoretical ansatz considered 
and the details of the Monte Carlo simulation used to 
generate the e^e" -^ Z^ samples. We then describe the 
proposed method, detector simulations and the imposed 
cuts for the event selection. Finally we present the fit 
techniques and the results obtained from Monte Carlo 
studies. 



Z$ SAMPLES 



Events of the signal e+e^ -^ ZH were generated us- 
ing the Pythia program [lj|. The cross section for the 
Higgsstrahlung process is given by: 



^ The new world average of the expected Higgs mass of 117 GeV 
Q is yet accesible in the current run of the Tevatron 
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with /3 = ^[s - (Mh + Mz)^][s - (Mh - MzY]/s. 

The effects of initial state brcmstrahhmg were included 
in the Pythia generation. 

For the production of a Higgs boson $ with arbitrary 
CV properties, e+e" -^ Z$, the amplitude for the pro- 
cess (which is not included in Pythia) can be described 
by adding a ZZA coupling with strength 77 to the SM 
matrix element as in ref.[l|. The matrix element of the 
process containing both the CV even amplitude, M.zh, 
and a CV odd amplitude, MzA, is given by: 
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FIG. 1: Definition of the production and decay angles of the 
process e+e" — > Z$[Z — > //]. 



where 77 is a dimensionlcss factor. The total cross section 
depends of the value of rj as follows [15| : 



following ratio: 
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and Mz, Fz denote the Z boson mass and width, Gp 
is the Fermi constant, and Ve^ CLe are the usual vector and 
axial-vector coupling constants of the e to the Z boson. 
In the SM 77 is zero. In the MSSM a ZZA is forbidden 
at Born level, but is induced via higher-order loop effects 
[9j. In general in extensions of the Higgs sector, 77 need 
not to be loop suppressed, and may be arbitrarily large. 
Hence, it is useful to allow for 77 to be a free parameter 
in the data analysis. 

As it was mentioned in the Introduction, the quantum 
numbers J^*^ of the Higgs bosons can be determined 
at future e^e~ linear colliders in a model independent 
way by analysing the angular dependence of the Hig- 
gstrahlung process. The most sensitive kinematic vari- 
able to distinguish the different contributions to Higgs 
boson production is 6*, the polar angle of the Z boson 
w.r.t. the beam axis in the laboratory frame. The sensi- 
tivity can be increased by including the angular distribu- 
tions of the decay to fermions, Z ^ ff in the boson rest 
frame (See Fig.^. Here the z-axis is chosen along the di- 
rection of the Z-boson momemtum. The decay amplitude 
is then a function of the angle between the Z momentum 
and /, 9* , and the angle between the Z production plane 
and the Z decay plane, (f>*. 

To obtain the angular distributions corresponding to 
the non-SM Higgs $ with arbitrary CV properties in 
the process e+e" -^ Z$, we have used a "re-weighting" 
method. This procedure allows one to obtain the dis- 
tributions for arbitrary values of 77 by weighting the dis- 
tributions for 7/ = according to the differential cross 
section in 0, 0*,0*. The weight factor is given by the 
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The squared amplitude |A^z$(?y)P has three contribu- 
tions: 

|Xz$(77)P - \Mzh\^+V^^MM*zhMza)W\Mza\^ 

(6) 
The first term reproduces the SM-like cross section. The 
interference term between the CV even and CV odd am- 
plitudes, linear in 77, generates a forward-backward asym- 
metry, that is a hallmark of CV violation. The third term 
correspond to the pseudoescalar Higgs cross section. Of 
course, 7; = brings us back to the scalar SM Higgs 
production. The explicit expression for the squared am- 
plitude of the brcmsstrahlung process in terms of 9, 9*, 
(j)* is taken from ref. 10] . In our procedure, the weight is 
re-scaled to be lower than 1, for a better treatment of er- 
rors. To check the reliability of the method we compared 
the obtained distributions using Monte Carlo with the 
analytical expressions. Figure |5]shows the obtained pro- 
duction angular distribution for the process e^e~ —> ZA 
using the procedure described above along with the ana- 
lytical form. The distribution is in very good agreement 
with the theoretical expectation proving the validity of 
the "reweighting" procedure. 



DESCRIPTION OF THE METHOD AND MONTE 
CARLO STUDIES 

We consider the production of Higgs events at the 
TESLA operating at a center-of-mass energy of 350 GeV, 
assuming an integrated luminosity of 500 fb~^. At this 
energy the main production process for the Higgs boson 
in the SM is the Higgsstrahlung process, e~^e~ — > ZH 
|l9l |. The corresponding expected number of events for 
this process is 6.6 x lO'^. 

We have choosen for the present study the process 
e+e" -^ ZH -^ fj,~^^~H with a Higgs boson mass of 




%Wf¥/^*' 




-0.2 0.; 
cos G 




115 116 



17 1ia 119 120 121 122 12J 124 125 



FIG. 2: Angular distribution of the process e e~ — > ZA ob- 
tained by the "re-weighting" method (See text) for an inte- 
grated luminosity of 500 fh~^ and a center of mass energy of 
^/s = 350 GeV, assuming a Higgs mass of 120 GeV. The line 
indicates the exact theoretical dependence. 



120 GeV. This decay channel exhibits a clean signature 
in the detector and the selection efficiencies are expected 
to be independent of the decay mode of the Higgs bo- 
son. We allow the produced Higgs to be either scalar or 
a mixture state $ including an interference term. 

All the Monte Carlo samples have been generated with 
the Pythia program as described in the previous sec- 
tion. These events are then passed through the simu- 
lation package SIMDET |2(]j, a parametric Monte Carlo 
program for a TESLA detector [21| which follows the pro- 

Bisal presented in the TESLA Conceptual Design Report 
. For the Higgs boson all decay modes are simulated as 
expected in the SM. The following background processes 
are considered in the analysis: e^e~ — » e^e^/^/^, 
e+e- -^ f+f~ij), e+e- -^ W+W' and e+e" -> ZZ. 
Both signal and background events are processed by the 
detector simulation package. 

For the event selection we follow ref . [15| . At least one 
muon and anti-muon are identified, with energy larger 
than 15 GeV. The mass of the di-muon system is re- 
quired to be consistent with the Z boson hypothesis 
within 5 GeV. The recoil mass of the di-muon system 
Ml, = (yi-(^M+ +E^-)f-iP^+ +P^- f has to be con- 
sistent with the H boson hypothesis within 5 GeV. This 
variable will yield a peak for the signal of the Higgs boson 
mass, independently of the Higgs boson decay mode. To 
remove a significant part of the remaining background, 
the absolute z-component of the di-muon system is re- 
quired to be smaller than 120 GeV. 

The momemtum of the selected muons are used to cal- 



FIG. 3: Recoil mass spectra off the Z in e e^ -^ ZH -^ 
fi'^ ^~X. Shadowed area represents the expected background 
events. 



culate the cosines of the production and decay angles for 
futher use in the method to determine the J^'-^ properties 
of the Higgs boson. It has been noted in []| that having 
excellent momemtum and energy resolution will allow the 
Z to be well reconstructed. The recoil mass against the Z, 
can then be used to detect the Higgs boson and to study 
its properties. Figure O shows the recoil mass distribu- 
tion for the rriH = 120 GeV signal, obtained from the 
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selected events in the sample of e 
The Higgs boson signal appears on top of a small back- 
ground. In Figure Qlthe corresponding cosO distribution 
is shown. The expected background is also presented. 
The combination of the cut on the z-component of the 
di-muon system, and the decreasing muon identification 
performance results in an efficiency for cosO > 0.9 close 
to zero. 

The kinematics of the e+e~ -^ Z^ [Z -^ //] process 
is described by the production and decay angles 0,9*, 
(j)* . The method we propose consists in generating 3- 
dimensional distributions in cos 6, cos 9* and cos cf>* using 
the Monte Carlo events generated as described above for 
each contribution in equation (3). We write then the 
likelihood: 



^= n 

(cos 6)i,(cos 0*)j .(cos 0*)fc 



f-^ijk 
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(7) 



where Ndata {hj^k) is the number of events of the hypo- 
thetical data sample and pujk is the expected number in 
the ijk-th bin. fiijk is calculated assuming a linear com- 
bination of the number of events of three Monte Carlo 
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TABLE I: Obtained values of the parameters for differents 
"data" samples 



FIG. 4: Angular distribution, cos 9, of the selected events in 
e^e~ -^ ZH — > fj,'^ fj,~X. The shaded histogram correspond 
to the expected background. 



samples, corresponding to the production of scalar Higgs 
(MC_ZH), pseudoscalar (MC_ZA) Higgs and events for 
the interference term (MC_IN): 

fi,,k = M.ia.MC.ZH,,k+[3-MCJN,,k+-f-MC^A,k) 

(8) 
where Af is the overall normalization factor between num- 
bers of data and Monte Carlo events which can be fixed 
{Af = 0.1 in our case) or left free as a further check of the 
fit. The likelihood is then maximized with respect to a, /3 
and 7. The absolute value of (3 indicates the contribution 
of interference term in the sample and a and 7 indicate 
the fraction of scalar and pseudoscalar components re- 
spectively. A significant deviation of P from zero would 
imply the existence of CV violation, independent of the 
specific model. For a scalar Higgs sample (77 — 0), the 
result of the fit is expected to be a = 1 and (3 — ■^ — 0. 

We have performed Monte Carlo studies with several 
hypothetical data samples with non-standard values of 77. 
A maximum likelihood fit for the best linear combination 
of MC.ZH, MC.ZA and MC.IN to match the hypothet- 
ical data sample gave statistical errors of 0.04, 0.02 and 
0.04 for a, /? and 7, respectively. The results of these 
studies using different values of rj are given in table H 

The value of a gives the fraction of the scalar J^'~' = 
0++ component of the Higgs boson, while 7 gives the 
contribution of the pseudoscalar Higgs component and 
increases quickly with 77 as expected. It can be seen from 
our results that the Monte Carlo study using a sample of 
pure scalar SM-like Higgs gives a consistent answer. This 
indicates the high sensitivity of the method to distin- 
guish a purely CT^-even state from a pseudoscalar CV-odd 
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-0.4 


0.002 ± 0.03 


-0.05 ± 0.02 


0.98 ± 0.04 


-0.25 


0.08 ± 0.04 


-0.06 ± 0.02 


0.92 ± 0.04 


-0.1 


0.43 ± 0.04 


-0.09 ± 0.02 


0.57 ± 0.04 


-0.05 


0.69 ± 0.04 
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0.97 ± 0.05 


0.003± 0.02 


0.03 ± 0.04 


0.05 


0.70 ± 0.05 


0.05 ± 0.02 


0.29 ± 0.04 


0.1 


0.40 ± 0.04 


0.04 ± 0.02 


0.59 ± 0.04 


0.25 


0.08 ± 0.04 


0.04 ± 0.02 


0.92 ± 0.04 


0.4 


0.002 ± 0.03 


0.01 ± 0.02 


0.98 ± 0.04 



state. Secondly, the method also allows one to determine 
whether the observed Higgs boson is a CV mixture and, 
if so, measure the odd and even component. It is evident 
that the statistical uncertainties prevent us to a large ex- 
tent from measuring the interference term. It should be 
noted that for Z — > fi^fi", as well as for Z — > e"*'e~, 
the interference term is suppressed by the smallness of 
Vf independently of the size of r/. However, the simulta- 
neous existence of fractions a and 7 would indicate CV 
violation for the ZZcjj coupling. The method proposed 
here gives sensible results in the case that there is any 
significant CT'-even component in the Higgs boson or if 
(j) is almost purely CV-odd. The statististical significance 
can certainly be increased including the e'^e" — > e'^e^ X 
channel. 



SUMMARY 

We have proposed a novel method for the measurement 
of the parity of the Higgs boson using the angular distri- 
butions of the differential cross section of e+e~ — > Z(f). 
The statistical power of our method using Monte Carlo 
generated hypothetical data samples is shown in Table 
^ The results indicate that, for an integrated luminosity 
of 500/5-1, g^^ 35Q QgY centre-of-mass energy, TESLA 
will be able to unambiguosly determine whether a Higgs 
boson is a state 0++ (CT'-even, scalar) or has a contri- 
bution of the ^ (CT'-odd, pseudoscalar) state, like in 
general extensions of Higgs model. We also estimate the 
statistical uncertainties for the measurement of the CV 
violating interference term. We hope that this technique 
will allow confirmation of the expected J^ assignment 
of a Higgs boson candidate. 
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